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SHOCKED SIMILARITY COLLAPSES AND FLOWS IN STAR FORMATION PROCESSES 
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ABSTRACT 

We propose self-similar shocked flow models for certain dynamical evolution phases of young stellar objects 
(YSOs), 'champagne flows' of H n regions surrounding OB stars and shaping processes of planetary nebulae 
(PNe). We analyze an isothermal fluid of spherical symmetry and construct families of similarity shocked flow 
solutions featured by: 1. either a core expansion with a finite central density or a core accretion at constant rate 
with a density scaling oc r~ 3 / 2 ; 2. a shock moving outward at a constant speed; 3. a preshock gas approaching 
a constant speed at large r with a density scaling oc r~ 2 . In addition to testing numerical codes, our models can 
accommodate diverse shocked flows with or without a core collapse or outflow and an envelope expansion or 
contraction. As an application, we introduce our model analysis to observations of Bok globule B335. 

Subject headings: Hn regions - HD - ISM: clouds - shock waves - stars: formation: B335 - winds, outflows 



1. INTRODUCTION 

Far away from initial and boundary conditions, dynamical 
evolution of a fluid may lead to self-similar phases with variable 
profiles shape-invariant and magnitudes properly scaled (Sedov 
1959). For star formation, self-similar solutions were found to 
describe collapses of isothermal gas clouds (Larson 1969; Pen- 
ston 1969; Shu 1977, hereafter S77). From an initial time t = 
when the core becomes singular to the final stage t — » oo, Shu 
(1977) derived the 'expansion-wave collapse solution' (EWCS) 
that has a free-fall core of density p oc r -3 / 2 , a radial infall speed 
// oc r" 1 / 2 , a constant mass accretion rate at r — > and a static 
envelope of a singular isothermal sphere (SIS) with the stagna- 
tion point moving outward at the sound speed a. This EWCS 
scenario of 'inside-out' collapses for star formation has been 
advocated by Shu et al. (1987) and compared with observa- 
tions (e.g., Zhou et al. 1993, hereafter Z93; Choi et al. 1995; 
Saito et al. 1999, hereafter Sa99; Harvey et al. 2001). 

Lou & Shen (2004, hereafter LS04) re-examined this prob- 
lem and derived new solutions in the 'semi-complete space' 
(0 < t < +oo), in contrast to the 'complete space' (-oo < t < 
+oo) of Hunter (1977) and Whitworth & Summers (1985, here- 
after WS). Gas flows at large r together with the two distinct 
asymptotic behaviors near the origin give rise to infinitely many 
similarity solutions. Those solutions of 'envelope expansion 
with core collapse' (EECC) are especially interesting. 

Besides these shock-free solutions, Tsai & Hsu (1995, here- 
after TH95) constructed a shocked self-similar collapse solu- 
tion for a situation where a central energy release initiates an 
outgoing shock during a protostellar collapse of a low-mass 
star. Such a shocked solution, matched with a static SIS en- 
velope, has a free-fall core yet with a lower mass accretion rate 
compared with the EWCS. TH95 also constructed a shocked 
expansion solution of a finite core density matched with a static 
SIS envelope, which was recently generalized and subsumed 



into a family of shocked 'champagne flow' solutions by Shu 
et al. (2002, hereafter S02) for expansions of H u regions 
around OB stars (Strtimgren 1939; Tenorio-Tagle 1979, 1982; 
Newman & Axford 1968; Mathews & O'Dell 1969; Franco, 
Tenorio-Tagle & Bodenheimer 1990, hereafter FTB). 

The main thrust of this Letter is to show that in the solution 
framework of LS04, a variety of shocked similarity flow solu- 
tions can be constructed and applied to astrophysical systems. 

2. SHOCKED SIMILARITY FLOW MODELS 

To study basic properties of shocked similarity flows of 
spherical symmetry, we assume isothermality with a constant 
sound speed a. The standard self-similar nonlinear ordinary 
differential equations (ODEs) (S77; TH95; S02; LS04) are: 
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m = x 2 a(x-v) , (3) 

where the independent similarity variable x = rj(at) and a(x), 
v(x), m(x) are the reduced density, radial speed and enclosed 
mass, respectively. The physical density p, radial speed u and 
enclosed mass M are obtained by the similarity transformation 
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u{r, t) = av(x) , M(r, t) = —m(x) . (4) 
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The analytical asymptotic solutions are (i) for x — > +oo, 

v — ► V , a — > A/x 2 , m—¥Ax, (5) 

where V and A are two parameters; and (ii) for x — ► 0, either 

v — ► — (Into/x) 1 ' 2 , a — > [mo/(2x 3 )]'/ 2 , m-^mo, (6) 
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where mo and B are two constant parameters (e.g., LS04). 

There exist two exact solutions, viz., the SIS solution (Bon- 
ner 1956; Chandrasekhar 1957) and the homogeneous 'Hubble 
flow' solution (WS; S02). ODEs (l)-(3) can be solved numer- 
ically to connect proper asymptotic solutions (5) - (7). 

The isothermal shock conditions 1 are the mass conservation 

a d (v d -x sd )a d = a u (v u -x su )a u 

and the momentum conservation 

a d [l + v d (v d - x sd )]a 2 d = a u [ 1 + v„( v u - x su )]a 2 u 

with the energy conservation involving radiative losses (e.g., 
Courant & Friedrichs 1976; Spitzer 1978), where subscript d 
(u) denotes the downstream (upstream) of a shock, a d x sd = a u x su 
is the shock speed, a d and a u are the downstream and upstream 
isothermal sound speeds, and x d = r/(a d t) and x„ = rj(a u t) are 
the downstream and upstream independent similarity variables. 
For a sound speed ratio r = a d j a u with rx sd = x su , we have 



v d -x sd -T(v u -x su ) = (TV d -v u )(v u -x su )(v d -x sd ) , 
a d /a u = (v u -x su )/[T(v d -x sd )] . 



(8) 



The special isothermal shock conditions of r = 1 with x sd = 
x S u = x s (TH95; S02) are adopted in our numerical examples. 

We take the so-called type 2 solutions (Hunter 1977, 1986; 
LS04) that cross the sonic critical line x— v = 1 once at x* < 1 
and approach asymptotic diverging behavior (6) for x — > as 
downstream (post-shock) flows. By integrating a type 2 solu- 
tion toward the sonic critical line, imposing shock condition (8) 
at each integration step and continuing further toward x — > +oo 
(upstream), we construct a family of shocked similarity flows 
displayed in Fig. 1 with r = 1. The resulting upstream solutions 
approach asymptotic behavior (5) as x — > +oo. Such a solution 
is uniquely determined by two parameters: the central mass ac- 
cretion rate mo (or equivalently, the crossing point x* < 1) and 
the shock location x s . Simultaneously determined are the veloc- 
ity and density parameters V and A as x — > +oo in (5). For each 
such solution (denoted as Class I), a core collapses with a con- 
stant mass accretion rate moa 3 /G and a shock travels at speed 
ax s outward to match with either an expanding or a contracting 
envelope; at very large r, the flow approaches a constant speed 
as either wind (including breeze and SIS) or inflow - all with 
density scalings of oc r~ 2 (LS04). 

For downstream solutions with asymptotic behavior (7), we 
can construct Class II solutions of shocked flows in parallel 
with the Class I solutions. Such a solution is uniquely deter- 
mined by the finite core parameter B and the shock location x s . 
For each Class II solutions, a core expands and a shock trav- 
els outward at a speed ax s to match with envelope expansion or 
contraction of a constant speed V ^0 or a breeze V = (S02) as 
x — > +oo shown by heavy solid curves in Fig. 1. 



Type 


Description 


Shock location x s 


V 


A 


Class I 


m = 0.406, z» = 0.23 


0.43 


-2.495 


0.250 






1.03 


-1.126 


0.847 






1.37 


0.0 


1.865 






1.43 


0.235 


2.158 


Class II 


8 = 0.1 


1.30 


-2.192 


0.0263 






2.00 


-0.711 


0.148 






2.80 


0.637 


0.640 



A few important features are noted here. 1. Shocked simi- 
larity solutions of Class I qualitatively share properties of the 
EWCS (inside-out collapse) at small r with density p oc r" 3 / 2 , 
infall speed u oc r~ l l 2 and constant core mass accretion rate mo. 
2. The main distinction is that the upstream solutions for both 
classes can be an inflow or a wind (LS04) or a breeze (S02) or 
a SIS 2 (TH95). 3. For very large B or small m (or small jc„), 
the downstream solution v(x) can oscillate (Hunter 1977 for the 
former and LS04 for the latter). The core mass accretion rate 
can be extremely low for the latter. 
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FIG. 1 . — Examples of Classes I and II (light and heavy solid curves) shocked 
flow solutions. For Class I solutions, we choose the downstream solution to 
cross the sonic critical line x-v = 1 [the dash-dotted line for v(x)] at x f = 0.23 
with mo = 0.406. An infinite number of solutions can be constructed for var- 
ious shock locations x s to match asymptotic behavior (5) as x — > +oo with 
various V and A (either a wind or a breeze or a contraction). This special 
downstream solution, referred to as the CP1 EECC solution in LS04, can an- 
alytically cross the sonic critical line again at x*(2) = 1.65 (dashed line). For 
Class II solutions in parallel, we choose the downstream solution of B = 0. 1 
and various x s to match (5) with various V and A. Parameters of these examples 
are summarized in Table 1. The dotted lines stand for the EWCS. 



3. STAR FORMATION PROCESSES 

3.1. Protostellar Collapses and Outflows ( Class I solutions) 

Regarding the EWCS scenario for star formation (Shu et al. 
1987), TH95 argued that a central energy output during a pro- 
tostellar collapse can initiate a shock into a SIS envelope with 
a lower core mass accretion rate. Based on solutions of LS04, 
we envision a much more general scenario for a central inside- 

1 The two jump conditions in the shock framework are the mass and momentum conservations pd(uj-u s ) = p u (u u -u s ) and a 2 d pd+ pdUj(u d -u s ) = a^p u + p u u u (u u -u s ), 
where u s is the shock speed. 

2 There are infinitely many discrete downstream solutions with extremely large B or small m Q and oscillating v(x) around zero to drive shocks into a SIS envelope. 
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out collapse with an adjustable central mass accretion rate and 
with an outgoing shock into various possible envelope flows 
such as winds, breezes, inflows or SIS at large r. The core 
mass accretion rate is M(0,f) = mocr'/G. While the EWCS has 
mo = 0.975 and the example of TH95 has mo = 0. 105, our Class 
I solutions for shocked flows give rise to a continuous range of 
< mo < 0.975. These solutions can describe certain similar- 
ity evolution phases during protostellar collapses and outflows. 
With these physical concepts, we now turn to the well-observed 
protostellar cloud B335 that possesses a collapsing core. 

Besides the observed bipolar CO outflows on large scales, 
a SIS environment is presumed to surround the central col- 
lapse region and has been modeled using the EWCS (S77; Z93; 
Sa99). From a temperature profile estimated from dust emis- 
sions (Zhou et al. 1990), B335 globule was approximated as 
an isothermal cloud with an effective sound speed a e ff ~ 0.23 
km s -1 . In essence, a gross spherical symmetry is assumed for 
B335 cloud by regarding the bipolar outflows as additional fea- 
tures produced by a circumstellar disk within < 100 AU (Har- 
vey et al. 2003). In the EWCS modeling (Z93), a core collapses 
as a free fall while an envelope remains static with an estimated 
infall radius r m f ~ 0.03 pc. 

Given some successes of the EWCS model in interpreting the 
density and velocity profiles (via simulations of spectral lines; 
Z93; Choi et al. 1995; Velusamy et al. 1995), departures from 
the EWCS model exist. In Sa99, the estimated infall speed at 
r ~ 2200 AU is - 0. 14 km s" 1 for r mf ~ 7200 AU, smaller than 
the EWCS prediction of ~ 0.34 km s" 1 by a factor of ~ 2.4, 
and the estimated mass accretion rate is 1.2-2.4 x 10~ 6 M Q 
yr" 1 , corresponding to mo ~ 0.418-0.836 that tends to be 
less than the EWCS prediction m = 0.975 (Z93). Meanwhile, 
the column density of the outer envelope at r = 10,000 AU is 
N ~ 6.3 x 10 21 cm" 2 , slightly larger than the EWCS prediction 
N — 5.4 x 10 21 cm" 2 (Sa99). This last difference is aggravated 
by the near-IR extinction study of Harvey et al. (2001) that 
gives an envelope density 3-5 times higher than EWCS pre- 
diction while the density p scales as r~ 2 . HCN spectral analysis 
of B335 seems to indicate a central core collapse with an enve- 
lope expansion distinctly different from bipolar CO outflows. 

The major advantage of our shocked similarity flow model 
is the capability of accommodating observations of B335 glob- 
ule while making testable predictions for envelope expansion in 
contrast to a SIS of the EWCS. As an example, we introduce the 
unshocked CP1 EECC solution (dashed curves in Fig. 1) that 
analytically passes the sonic critical line twice (LS04). While 
the inner core freely falls with a p oc r~ 3//2 , the envelope expands 
with a p oc r~ 2 . In contrast to the EWCS, the transition across 
r m f is smooth without a 'kink' . Still taking x = 1 for r m f = 7200 
AU (or an age of the similarity process is ~ 1.5 x 10 5 yr in the 
EWCS model, Z93), the infall speed at r ~ 2200 AU is - 0.12 
km s" 1 and the mass accretion rate given by CP1 EECC solution 
is m = 0.406 less than a half of 0.975 in the EWCS. At large 
r, the density approaches p <~ Aa 2 ^/ (AirGr 2 ) with A = 5.158 

3 In reality, a polytropic model should be more pertinent as discussed in LS04 

4 Practically, it might be far easier to detect the presence of large-scale shocks 
signatures amenable to observations, i.e., through emission diagnostics from hij 



(Table 3 of LS04). That is, our CP1 EECC model gives an 
envelope density <~ 2.5 times higher than that of a SIS in the 
EWCS (i.e. A = 2). This tends to agree with the inference of 
3-5 times by Harvey et al. (2001). One distinctive conse- 
quence of our isothermal model is that the envelope, instead of 
a SIS, expands with a speed 3 of ~ 0.25 km s" 1 at r = 25,000 
AU. This is qualitatively consistent with the preliminary HCN 
spectral analysis (A. Lapinov 2004, private communications). 
By inserting a shock at a proper location, it has a range to ad- 
just the data fit for density and speed profiles and there is a 
preliminary evidence of shock in B335 system 4 (e.g., Nisini et 
al. 1999). Sharing the qualitative features of the EWCS in the 
central region, our Class I solutions (including CP1 EECC so- 
lution) of shocked similarity flows are versatile to model other 
protostellar cloud or Bok globule systems. 




FIG. 2. — The family of Class II solutions with an invariant downstream 
(S02) as B -> (B = 10 -6 ). Here, a(x) is rescaled as R(x) = x 2 a(x)/B. The 
heavy solid line stands for the 'champagne breeze' of S02. In addition to 
breezes and contractions at large r, many shocked similarity flows approach 
constant wind speeds as x — > oo that can be either supersonic or subsonic. 



3.2. 'Champagne Flows ' in Hu Regions ( Class II solutions) 

As massive OB stars form in molecular clouds, intense ul- 
traviolet radiations ionize surrounding Hi clouds and carve out 
Hu regions (Stromgren 1939). As an ionization front passes 
by, a pressure gradient develops between H n and H I regions 
to drive flows. Shocks can naturally occur in such an impul- 
sive expansion phase. An expansion of a Hu region embedded 
in a uniform ambient medium has been studied (e.g., Newman 
& Axford 1968; Mathews & O'Dell 1969). Other possible en- 
virons may involve power-law cloud density profiles of ex f~" 
with 1 < n < 3 and a mean of n = 2 (e.g. Arquilla & Gold- 
smith 1985). As ionization or shock fronts encounter a strong 
negative density gradient, the ensuing supersonic expansion of 



For a polytropic gas, the asymptotic wind speed should tend to vanish at large r. 
than to estimate the velocity and density profiles as shocks will leave non-kinematic 
;hly compressed post-shock gas or chemical or radiative cooling processes. 
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the ionized gas is referred to as the 'champagne phase' 5 (FTB). 
When a fast moving ionization front cannot be confined under 
certain conditions (FTB) and ionizes the entire cloud in a rela- 
tively short time, the inevitable expansion of Hn regions might 
evolve into a self-similar phase (S02). 

In the scenario above, a cloud is presumed static initially and 
the H ii region forms instantaneously at / = when the stellar 
nuclear burning turns on. Our more general scenario would al- 
low a molecular cloud with asymptotic flows (either winds and 
breezes or accretions and contractions) at large r. To model ex- 
pansions of Hn regions as the ionization front travels, we focus 
on various possible large-scale systematic flows at large r and 
put aside messy central star-forming processes on much smaller 
scales. Within the ionized cloud sphere, asymptotic winds and 
accretions at large r can be connected to inner champagne flows 
via transonic shocks. As transients peter out, an expanding H 
ii region may become self-similar with an outgoing shock to 
match with an asymptotic wind, breeze or accretion. 

The shocked breeze solutions has been used to model simi- 
larity expansions of Hn regions (S02). For each 'champagne 
breeze' solution of S02, we can use the same downstream to 
construct Class II solutions with shocked upstream flows which 
approach breezes, winds and accretions at large r. As B — > 0, 
the downstream solution becomes invariant (S02) and the re- 
sulting family of upstream solutions is displayed in Fig. 2. 

Observations do suggest the presence of systematic flows in 
Hn regions and models involving stellar winds were proposed 
(e.g. Dyson 1977; Hippelein & Munch 1981; Comeron 1997). 
Such outflows may be supersonic of a few km s -1 (e.g., Re- 
lano et al. 2003). In reference to our Class II models, shocked 
similarity 'champagne winds' in H n regions distinctly differ 
from 'champagne breezes' or SIS by having both constant wind 
speeds at large r and outgoing shocks at constant speed. Com- 



pared with breezes, 'champagne winds' would dynamically im- 
pact the surrounding interstellar medium more efficiently. 

3.3. Implications for Planetary Nebulae ( PNe ) 

In hydrodynamic models for shaping PNe, the interacting 
stellar wind (ISW) scenario (Kwok et al. 1978) is often in- 
voked. In a few spherical cases, a much faster stellar wind from 
the core catches up with a slower dense wind - the remnant of 
the AGB phase, and form shocks. Our shocked similarity flow 
solutions of Class I clearly show the physical feasibility that the 
shaping of a PNe involving ISW can be concurrently accompa- 
nied by a central accretion or infall towards a proto white dwarf 
in a similarity manner. In order to model bipolar PNe other 
than round or mildly elliptical morphologies, aspherical aspects 
are included in the so-called generalized ISW (GISW) model 
(Balick & Frank 2002). By flow instabilities, a shocked enve- 
lope can expand to create various nebula morphologies while 
the central core continuously accretes materials to form a white 
dwarf. As a proto white dwarf continuously accretes materials 
to approach the Chandrasekhar limit of ~ 1.39M , there might 
be even a possibility of igniting a type la supernova explosion. 
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